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= ABSTRACT 

Experience gained from the dynamic analysis of horizontal axis turbine rotors based on the use of 
the HCSTAB-HFW computer code is described. Three topics are covered, dealing with the frequencies of a 
— rotating beam, the use of the fundamental mode of a uniform cantilever beam, and the analysis of reso- 

nance dwell. Immensely high peak loads were generated by the code for resonance dwell , indicating 

r further need for including structural damping and for transient analysis capability. The effect of 

structural damping, newly incorporated In the code, is finally described.-. 


INTRODUCTION 


The ^K)STAB-HFW computer code was developed 
by Paragon Pacific under support from NASA Lewis 
Res^rch Center for the dynamic analysis of 
ho r 1 zon t a 1-axis wi nd turb ine ro tor s . Expe r 1 en c e 
on the use of the code is described In this paper 
Three topics are covered. 


FREQUENCIES OF ROTATING BEAM 

rihe first topic deals with ..the frequencies 
of a rotating beam, which are needed as input to 
the MOSTAB-HFW code. Based on the numerical ^ 
results of Yntema (Ref. 1), a simple expression 
has been derived relating the frequency of a 
rotating beam to that "dT the same beam without 
rotation. For wind turbine blades with commonly 
used dfmen^h s , t he ex press i bn is very co nven - 
lent and accurate enough for preliminary design 
and analysis purposes. The derivation is pre- 
sented below. 


Yntema (Ref. 1) has given the frequency of 
a rotating beam by the equation 

where and frequencies of the 

beam with and without rotation, respectively, 
and are dimensionless coefficients, the 
subscript n~ denotes the mode, e is the ratio 
between the offset and the beam length, and 
the rotation speed. This equation is quite gene- 
ral and applicable to nonuniform beams with a 
large variety of mass and stiffness distributions. 
For linear mass and stiffness distributions 
Yntema calculated Kon and Kin the first 
three modes of hinged and cantilever beams. The 
linear distributions vary from a beam with 
constant-mass and/or constant-stiffness to a 
beam whose mass and/or stiffness diminishes all 
the way down to zero at one end of the beam. 

Yntema ‘s results for Kqi and for the 

fundamental modes of cantilever beams deserve 
special attention, because these modes play an 
important role In the dynamics of wind turbine 
blades and are needed 1n the MOSTAB-HFl^ code. 

Such results are given In his Figs. 15 and 16. 


It is not worthy that, for the wide range of 
linear mass and stiffness distributions these 
coefficients lie within limited ranges: 

1,17<Kq^< 1.26, 1.57<Kj^<1.93 

Thus, Kqi already lies within a very narrow band. 
Although Kn does not, we note that the offset 
ratio e is seldom greater than 0.1 for wind 
turbines. Hence, by assuming a maximum value of 
e " 0.1, we have 

0.157 < < 0.193 

and, for e<0.1 

1.17<(Koi + K^^e) < 1.45 

This interesting result shows that (Kqi + ^11^) 
has only a limited variation for e<0.1. Within 
11 per cent accuracy w^“prDpose to take simply 

(Kqi + Kj^e) = 1.31 for e<0.1 

The fundamental frequency of a rotating 
cantilever beam Is then given by 

^Rl^ " ^NRl^ + 1.31 a 

For some blades analyzed by the use of this 
equation, the accuracy on has been found 
to be within 1 per cent. 

USE OF FUNDAMENTAL 
CANTILEVER MODE 

The next topic is associated with the natural 
modes of vibration of the blade which are also 
needed as input to the MOSTAB-HFW code. Although 
the code can accommodate coupled modes including 
flapwise and chordwise bending together with 
twisting of the blade, the use of the fundamental 
mode of a uniform cantilever beam for both flap- 
wise and chordwise bending has been found to be 
satisfactory for quick calculation of the dynamic 
loads. This 1s essentially to apply the classical 
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Rayleign procedure 1n structural dynamics 
to the computer computation. 

The use of the fundamental can ci lever mode 
simplifies the execution of the MOSTAB-HFW code 
greatly. The mode Is well known and can be 
entered as input once for all, regardless of 
the values of the associated bending frequencies. 
The advantage of this is obvious, particularly 
during preliminary design and analysis of the 
rotor blade, when many changes may be needed. 

The use of the fundamental cantilever mode 
in the code has been compared with the use of 
coupled modes in sample calculations. Differ- 
ences between the resulting blade peak and cyclic 
loads ranged from 8 per cent for the twisting 
moment to 23 per cent for the flapwise bending 
moment. These are considered satisfactory, 
since the loads generated by the MOSTAB-HFW 
code can deviate from the test data by as much 
as 25 per cent, as reported by Spera (Ref. 2). 


RESONANCE DWELL 

The third topic in this paper deals with 
resonance dwell. This is when the ratio between 
a natural frequency of the blade and the rotor 
speed is an integer. Such a situation is natur- 
ally to be avoided in the design for the normal 
operation of a wind turbine. However, it cannot 
be avoided during a transient operation such as 
start-up, shut-down, or emergency shut-down, when 
the rotor speed varies and passes through reso- 
nances with the blade. Resonance dwell represents 
the limiting case of constant rotor speed or zero 
rotor acceleration and can be handled by the 
MOSTAB-HFW code, although the code does not treat 
the general transient case with a non-zero 
acceleration. 

The importance of accommodating transient 
operations in the design of wind turbines cannot 
be overemphasized. Dynamic loads on the blades 
of the MOD-OA wind turbines were measured during 
actual transient operations and found to be four 
to five times the loads during normal operation 
(Ref. 3). High transient loads have also been 
reported by Dugundji (Ref. 4), who tested a small 
wind turbine model in a wind tunnel at MIT. 

It had been hoped that, in the dynamic 
analysis of a wind turbine, results of trim 
analysis during resonance dwell based on the use 
of the MOSTAB-HFW code would shed some light on 
what might happen when the rotor speed was 
varying and passing through a resonance. Unfor- 
tunately, the peak loads generated by the code 
for resonance dwell turned out to be too high 
to be realistic. For instance, in one case of 
resonance dwell the code generated a maximum 
flapwise bending moment in the blade 36 per cent 
higher than that for normal operating speed, a 
maximum chordwise bending moment more than 5 
times higher, and a maximum shaft torque more 
than 16 times higher. 

The very high loads for resonance dwell 
given by the MOSTAB-HFW code were due to lack 
of structural damping in the code. The results 
reflect the extreme sensitivity of the structurally 


undamped dynamic system as treated by the original 
code. They point to the need for including struc- 
tural damping and for transient analysis capabi- 
lity. 

EFFECT OF STRUCTURAL DAMPING 

In view of its apparent importance, struc- 
tural damping was recently incorporated in the 
MOSTAB-HFW code by Paragon Pacific, in addition 
to the already present aerodynamic damping. 

The same case of resonance dwell mentioned before 
was analyzed again by means of the code, only 
now with 2 per cent of critical damping added 
for both the flapwise and chordwise bending 
modes. The maximum flapwise and chordwise bend- 
ing moments in the blade and maximum shaft torque 
were found to be 69, 18, and 12 per cent, respec- 
tively, of their previous level for zero struc- 
tural damping. The load reduction was, therefore, 
drastic. 
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QUESTIONS AND ANSWERS 


Y.Y. Yu 


From: G* Beaulieu 


Q* your constant K values consider the reducing stiffness of the cantilever blade? 

If not , your resonant frequencies will be shifted, 

- A: The values of and given by Yntema aover a large variety of linear stiffness 

(as well as mass) distributions. The final approximate frequency equation in the 
paper therefore applies to cantilever blades with such stiffness variations , 

Fromr T, Currin 

Q: After rough approximations of input, is using a large code, i.e., MOSTAB, justified? 

A: The use of approximate frequencies and mode shapes as indicated in this paper saves 

time for preparing input during the early stage of design when a number of configur- 
ations may be considered and many changes needed. When the design is narrowed down^ 
such approximations are replaced by more accurate values while many other inputs to 
'^b.e code remain the same. The use of the code at an early stage therefore does not 
add much to the analysis work. 


From: Jack Landgrebe 


Q: Does your sensitivity demonstration indicate that you can get any answer with dif- 

ferent dampings and what answers do you believe? 

A: L^<ata on structural damping of wind turbine blades are lacking because little atten- 

^ paid to such damping. When it is better known, the transient loads 

aan be estimated more accurately , 


From: " W.C, Walton 


Q: 


1) Should results be compared with "normal operating loads" or with design loads? 

2) Please explain the term fixed axis. 


A: 2) The transient loads are compared with normal operating loads because both must 

be taken care of as far as structural strength is concerned, although their 
fatigue effects are different due to the different numbers of cycles during 
the life of the wind turbine. 


2) The MOSTAB-fiFW code assumes a rotor with a fixed axis in space. In other 

... words, the rotor is assumed to be uncoupled from the other components of 

the wind turbine system. In contrast , the MOSTAS code considers the inter- 
action between the rotor and the other components. 

-From: Anonymous 


Q: How was damping incorporated into MOSTAB and what is required for input data to 

take this into account? 

A: Modal damping was introduced for each mode in the form of a viscous damping term. 

Damping ratio (percentage of critical damping) is assigned to each mode as input. 



From: Jack Hoffman 

Q: Why did the flap loads go up with damping? 

A: Thq flapwise bending moment did not go up. With 2 percent damping, the moment 

was 69 percent of its value for zero damping. 
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